Specific heat of single crystal MgB 2 : 
a two-band superconductor with two different anisotropics 
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Heat-capacity measurements of a 39 /ig MgB2 single crystal in fields up to 14 T and below 3 K 
allow the determination of the low-temperature linear term of the specific heat, its field dependence 
and its anisotropy. Our results are compatible with two-band superconductivity, the band carrying 
the smaller gap being isotropic, that carrying the larger gap having an anisotropy of ~ 5. Three 
different upper critical fields are thus needed to describe the superconducting state of MgE>2 . 

PACS numbers: 74.70.Ad, 74.25.Bt, 74.60.Ec, 74.25.Jb. 
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Shortly after the discovery of 40-K superconductivity 
in MgB 2 Q], its nature was intensively studied. The iso- 
tope effect soon indicated phonon- mediated pairing [Q. 
However, MgB2 cannot be considered as a classic conven- 
tional superconductor; indeed several studies point to the 
existence of a gap much smaller than the expected BCS 
value: band-structure calculations [^,^1, STM || spe- 
cific heat (C) @, |, |, [HJ penetration depth and 
various spectroscopic experiments 

The so-called "two-band" model [pTI , which considers 
two different s-wave superconducting gaps (a larger and 
a smaller gap, A CT and A ff ) on separate sheets of the 
Fermi surface (the a- and 7r-band, respectively), now 
seems widely accepted and explains most experimental 
results jl9|, |2^] , with few exceptions [^lj . The theory of 
two-band superconductivity was suggested soon after the 
development of the BCS theory p2[ , and some s-d met- 
als have shown hints of this phenomenon, generally in the 
form of small deviations from BCS predictions, e.g. in C 
experiments [|4|. In contrast, MgB 2 shows such clear sig- 
natures that it might serve as a textbook example. For 
instance, the electronic specific heat of MgB2 shows an 
excess by an order of magnitude at ~ T c /5, and presents 
an exponential dependence at low temperature (T) in- 
dicating the existence of a small gap on a large fraction 
(~ 50%) of the Fermi surface g |, g, M ffl]. At low 
T the effect of a small magnetic field (H) on the coeffi- 
cient of the linear C-term, ^(H) = lim^^o C*(iJ)/T, is 
dramatic. Classically, ^(H), a thermodynamic probe of 
the electronic density of states at low-energy, is expected 
to be linear in H up to H C 2 p3fl . However, for MgB2 an 
extreme non-linearity was observed, which was assigned 
to the existence of an additional, smaller gap |g|, ^, [To| . 

The only specific heat experiment on single crystal re- 
ported so far was focused on the superconducting tran- 
sition [^5[ ^6). Such experiments are delicate since the 
mass of MgB 2 crystals is usually smaller than 100 /ig, 
but are required to study the anisotropy. Previous stud- 
ies showed that the anisotropy of superconducting prop- 
erties depends on both field and temperature J27], B8[ . 
This puzzling behavior motivated the present study. 

We measured the low-temperature specific heat of a 
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FIG. 1: Typical curves of total specific heat, including that of 
addenda, for different magnetic fields (from bottom to top); 
a) 0, 0.1, 0.2, 0.5, 1, 2, 3, 4, 6, and 10 T along the c-axis; b) 
0, 0.1, 0.2, 0.5, 1, 3, 6, 8, 10, 12, and 14 T in the (ab)-planes. 
The straight lines are linear fits with a common slope. 



MgB2 single crystal with the field parallel and perpen- 
dicular to the boron planes. The results clearly show 
that the extreme non-linearity of ^{H), first observed 
in polycrystals j7|, [|, ^|, |To|] , is an intrinsic property of 
MgB2, related to the existence of two superconducting 
gaps. Moreover, these measurements reveal a dramatic 
variation of the effective anisotropy with the field, from 
1 to ~ 5. We interpret these results as a manifestation of 
the two-band nature of superconductivity in MgB2 , plus 
the dimensionality of the 7r-band (isotropic) and cr-band 
(anisotropic). 

The single crystal of MgB2 is similar to those described 
in Ref. |29[ . Its mass, 39 /ig, is small for a quantitative 
study of specific heat. For this purpose, we developed 
a miniaturized version of the relaxation calorimeter de- 
scribed in Ref. 0. A commercial Cernox bare chip |3^] 
was thinned down to 1 mg mass and used simultaneously 
as sample platform, heater and thermometer. Two 10 /im 
phosphor-bronze wires soldered onto the chip acted as 
mechanical support, thermal links and electrical leads. 
Only low-T studies were possible, since the heat capacity 
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of the platform, essentially of phononic origin (sapphire), 
increases rapidly with T. We thus restrained our study 
to the range 2 to 3 K. Previous experiments have shown 
that this range of temperature lies sufficiently above the 
domain of possible parasitic magnetic contributions, but 
low enough (~ T c /15) to allow reliable determination of 
7(i?) (see e.g. Ref. Q). Owing to the small value of 
C at these temperatures, the relaxation time was short 
(0.1-0.5 s), necessitating aperture times shorter than one 
power line cycle for the digitization of voltages. In order 
to suppress the 50 Hz noise, the signal was averaged over 
50 to 200 relaxations. The field dependence of the ad- 
denda represents at most 13% of the normal electronic 
linear term of the MgB2 crystal, and was calibrated by 
measuring a 196 /xg Ag reference sample. 

Figure [l] shows the specific heat for both orientations. 
It is impossible to compensate exactly for the small but 
unknown quantity of GE- varnish used to attach the sam- 
ple. Instead, in order to calculate ~f{H), we used the zero- 
field curve as a reference. It was shown on polycrystals 
that the electronic specific heat of MgB2 is negligible in 
zero field in this temperature range |8|, g, |I(J : the parallel 
shift of the C/T curves in Fig. [I] is thus a direct measure- 
ment of ~/{H). This shift was precisly determined by fit- 
ting each curve with a linear electronic term "f(H)T plus 
a f3T 3 lattice term, imposing the same (3 coefficient for 
all fields. The result for j(H) is plotted in Fig. ||for both 
orientations. With the field perpendicular to the boron 
planes, ^(H) saturates above 3-4 T, indicating that the 
sample has entered the normal state. This value is con- 

(c) 

sistent with previous reports giving fXpHw w 3 T for 
the upper critical field along the c-axis |2l|, ^(| |27], [sT| . 
This also determines the normal electronic linear term 
7„ = 0.76 ± 0.03 mJK _2 gat _1 , comparable with litera- 
ture values 0, ||, ^, |l0| . The normal state is not fully es- 
tablished when the field is parallel to the (ab)-planes, in- 
is higher than the maximum applied 

{ c f ] « 18- 
Therefore, 



dicating that 

field of 14 T. An extrapolation suggests uqH, 
22 T, in agreement with reported values [p2| 
the anisotropy of the upper critical field at T <C T c is 
established in bulk. 

The anisotropy of H C 2 in MgB2 has long been de- 
bated p2| , but was finally assessed by experiments us- 
ing clean sing le crystals §|, §(| [27], H Tw0 kinds 
of anisotropy can be considered: the anisotropy of the 
superconducting gap, and that of the Fermi surface. To 
model the former, some fc-dependence of the gap has to 
be assumed. Such an approach was attempted using a 
one-band model for MgB2 |3j|. The two-band model 
is an extreme case of such an anisotropy: depending 
on which Fermi sheet is considered, the superconduct- 
ing gap may be either large or small. The effect of an 
anisotropic Fermi surface on superconducting properties 
can be modeled by renormalizing the isotropic results by 
a function of the ratio of the effective masses, a method 
that was successfully applied to cuprate superconduc- 
tors [53 . Angst et al. used this approach for MgB2 p7| , 
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FIG. 2: Coefficient of the electronic linear term versus mag- 
netic field applied parallel (□) or perpendicular (•) to the 
boron planes. The long-dashed line represents the normal 
state contribution. The short-dashed lines are guides for the 
eyes. Inset: expanded view of the low-field region; here the 
long-dashed line represents the partial normal-state contribu- 
tion of the small-gapped band (see text and Fig. W) . 



but showed that the shape of the torque curves did not 
follow this simple one-band, anisotropic model over the 
whole range of field orientations. 

From Fig. ^, it is clear that it is not possible to normal- 
ize the effect of the field by using a constant anisotropy 
factor. For H < 0.5 T, the values of j(H) are al- 
most undistinguishable, irrespective of the field direction, 
showing the absence of anisotropy (see inset of Fig. ||). 
However, a field scaling factor of ~ 5 is needed to super- 
pose the j(H) curves just below saturation, i.e. near 
H C 2- An effective anisotropy T e ff can be defined as 
the scaling factor by which the field along the c-axis 
must be multiplied in order to merge both j(H) curves. 
Figure |^ presents T e ff versus H//(ab). At low field, 
T e tf « 1, since the effect of H does not measurably de- 
pend on the orientation; T e ff rapidly increases with H, 
and tends toward the anisotropy of H c2 , ~ 5. Such an 
£f-dependent anisotropy was also deduced from torque 
measurements 27 . The present analysis shows that the 
change of ^ e ff with H (and possibly with T) is an in- 
trinsic property of MgB2, having a bulk thermodynamic 
signature, even though the amplitude of this variation 
may depend on the physical quantity by which it is de- 
termined. We will argue in the following that, to explain 
specific heat measurements, both kinds of anisotropy are 
required: the band dependence of the gap, plus a renor- 
malization by the effective masses on the cr-band. 

Let us first focus on the shape of the ~f(H) curves. Sim- 
ilarly to previous data on polycrystals, "/(H) increases 
sharply at low fields. The present data show that this ef- 
fect does not depend on the orientation of the field, and 
point to the existence of a large fraction of normal elec- 
trons in the mixed state already at very low H. In a typ- 
ical s-wave, dirty type-II superconductor, j(H) mainly 
originates from the vortex cores and should be approxi- 
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FIG. 3: Effective anisotropy, defined as the ratio of magnetic 
fields applied in the (afe)-planes and along the c-axis yielding 
the same 7-value. Note that the choice of the abscissa is 
arbitrary: we chose H//(ab), but we could just as well plot it 
versus H//c or versus 'y(H). 




FIG. 4: Separation of both contributions to f(H): ^"(H) 
(straight lines) and ^(H) for the field applied along the c- 
axis (•) and in the (afc)-planes (□). The dashed line is a guide 
to the eyes for ^(H). 



mately proportional to H (i.e. the number of vortices). 
The strong non-linearity of ^(H) for MgB2 is atypical: 
for fioH — 1 T, approximatively half of 7„ is recovered, 
irrespective of the field direction. The effect is particu- 
larly striking for H//(ab): about 5% of H^ b) drives ~ 50% 
of the electrons into the normal state. 

The peculiar shape of 7(H) can be understood by 
considering the meaning of a vortex within a two-band 
model. Both bands supply the carriers and currents that 
define one vortex, and both bands contribute to the lo- 
cal density-of-states (LDOS) in one vortex core. Sim- 
ulations |35| and STM studies Q suggest that a core 
results from the superposition of two peaks in the LDOS 
with different diameters defined as usual by the coher- 
ence length £0- or £ v associated with each band. The 
sharper peaks with diameter £ CT lead to a conventional 
core contribution for the er-band, and overlap at H c i- In 
contrast, the large coherence length of the 7r-band gives 
rise to broad LDOS peaks, with a diameter much larger 
than ($o/27r_ff C 2) 5 . These giant cores, which have been 
observed by STM |J , contribute heavily to the total low- 
T DOS, as measured by j(H), until they start to overlap 
much below H c2 . At this point, the contribution of the 
7r-band to j(H) saturates. At higher fields, superconduc- 
tivity in the 7r-band is maintained up to H C 2 by coupling 
to the er-band, but the variation of "f{H) now comes from 
the cr-band contribution. The 7 (if) curves of Fig. || com- 
pare qualitatively well with these predictions. Note that 
the analysis of previous measurements on polycrystalline 
samples j?], ||, [l(| was complicated by angular aver- 
aging of 7(iJ//c) and ^(H/fab), leading to a flattening 
above 8 T. 

Within the present scheme, the data of Fig. |2] can be 
used to determine the anisotropy and T a of each sub- 
system separately. The abnormally fast initial increase 
of j(H) at low H is attributed to the giant vortex cores 
associated with 7r-band carriers, and, as shown by the 
inset of Fig. EL is isotropic; therefore — 1. On the 



other hand, the er-band contribution is responsible for the 
anisotropy of H C 2 , which we determined when describing 
Fig. 2; therefore r CT ps 6. This is consistent with the 
3D and quasi-2D character of the n- and cr-band, respec- 
tively, found in band-structure calculations ||, [ll| p0[ , 
and with the experimental determination by torque mag- 
netometry at low T . The effective anisotropy shown 
in Fig. |^ is a weighted average which depends on the rel- 
ative contribution of each band to the condensate for a 
given H and T: it varies smoothly from at low fields 
to r CT at high fields. 

Besides the anisotropy of each subsystem, we can at- 
tempt to determine the full shape of the unusual con- 
tribution of the 7r-band, *y w (H). For this purpose, we 
assume that the cr-band contribution, j a (H), follows a 
classical behavior p3| , increasing linearly from to 7^ as 
the field goes from to H^oi H^ b \ depending on the 
field orientation. The linear sections above 1 T in Fig. || 
correspond to this contribution of the cr-band; their ex- 
trapolation to H = points to the normal-state value 
7^ of the 7r-band contribution (7^ + 7^ = j n ). By sub- 
tracting 7 CT (i?) from ■y(H), we isolate 7 7r (ff) (Fig. §). 

The parameters H^f, H^ b \ 7^, and 7^ are determined 
from Fig. ||. By requiring that the saturation value of 
7 7r (-ff) is isotropic, we find that the parameters cannot 
depart by more than ±15% from the average fitted values 

pit" 1 

The j^/ln ratio « 55/45 « 
1, is consistent with independent fits of the zero-field 
specific heat |11 and penetration depth and with 
band-structure calculations [Q, [l9], |20|. Furthermore, 
7 7r (iJ) ap pea rs to follow quantitatively numerical simu- 
lations |(| |35| ■ This shows that a model consisting of two 
superfluids of nearly equal weight, having different gaps 
and different anisotropy, gives a consistent description of 
various properties of MgB 2 . 

The recent measurements of the thermal conductivity 
of single crystals by Sologubenko et al. fjlf were inter- 
preted in a similar way. At low fields, the heat conduc- 



H H£ = 3.5 T, Aio#c2 b) = 19 T, 7 £ = 0.42 mJK" 2 
and 7£ = 0.34 mJK^gaf- 1 . 
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tivity increases sharply, for all field directions. As dis- 
cussed for 'y(H), this is attributed to the contribution of 
normal 3D 7r-band carriers in the anomalously large vor- 
tex cores. At higher field, this effect saturates and the 
anisotropy rises, due to the increasing contribution of the 
quasi 2D tr-band. A similar approach may explain not 
only why the anisotropy obtained from magnetic torque 
experiment changes with H and T, but also why a single- 
anisotropy model does not fit the angular dependence of 
H c2 |2l], g8|. We emphasize, however, that this view is 
simplified and phcnomenological. A more elaborate de- 
scription for the interaction between the two bands was 
proposed p6, 37 for the anisotropy of penetration depth 



and its temperature dependence (which can be under- 
stood as a consequence of the rapid thermal depletion of 
the isotropic 7r-band condensate). 

The phenomenological two-band model of Ref. Jll[] 
also uses two independent superfluids to explain the 
zero-field C of MgB 2 . Inter-band coupling was intro- 
duced by considering that the "virtual" T c of the 7r-band 
(A w /1.76k B ~ 13 K) is brought up to that of the er-band; 
the peculiar features in C then originate from this super- 
conductivity "above T c " . Here a somewhat similar situa- 
tion arises: a "virtual" upper critical field can be defined 



for the 7r-band, H* 2 ~ (^oA^ )/(H 2 v 2 F ) |38 ; superconduc 



tivity persists up to H^or H^ 2 "', depending on the field 
orientation, but above H£ 2 the overlap of huge vortex 
cores drives the majority of the 7r-band electrons normal 
(superconductivity "above H c2 "). Three characteristic 
fields are thus needed to describe the physics of MgB2: 



(ab) 



the two upper-critical fields TJ^and i^2 b ' )assoc i a tecl 
with the anisotropic er-band, and a cross-over field H£ 2 
associated with the 7r-band. We estimate the latter one 
by linearly extrapolating the low field 'y(H) to the partial 
normal-state 7 for this band (~ 0.557„), by analogy with 
classic superconductors. The inset of Fig. || shows this 
construction: we find noH^ 2 ~ 0.3-0.4 T. These fields 
determine the values of the three coherence lengths of 



the system: for the anisotropic er-band £, 



(ah) 



10 nm 



(c) 

and Q- ~ 2 nm; for the isotropic 7r-band £ w ~ 30 nm. 
The latter value compares favorably with a direct mea- 
surement by STM giving 50 nm for the vortex-core di- 
ameter in the 7r-band ||. All three fields are needed to 
give a phenomenological description of the superconduct- 
ing properties of MgB2. Note that the ratio between 



H 



(ab) 
c2 



(c) 

and H c2 is driven by T a , i.e. the anisotropy 
of the effective masses of the cr-band, whereas the ra- 
tio between H { c2 } ~ (^ A 2 a )/(h 2 v 2 F ) and H£ 2 is given by 
~ A 2 /A 2 , since the Fermi velocities are comparable on 

both bands. The experimental ratio H^ 2 /H£ 2 ~ 10 re- 
ported here agrees with various estimations of A a / A n ~ 
3-3.5 § H |! 0, H H 
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